MULTI-MATRIX VECTOR COHERENT STATES 
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Abstract. A class of vector coherent states is derived with multiple of ma- 
trices as vectors in a Hilbert space, where the Hilbert space is taken to be the 
tensor product of several other Hilbert spaces. As examples vector coherent 
states with multiple of quaternions and octonions are given. The resulting gen- 
eralized oscillator algebra is briefly discussed. Further, vector coherent states 
for a tensored Hamiltonian system are obtained by the same method. As 
particular cases, coherent states are obtained for tensored Jaynes-Cummings 
type Hamiltonians and for a two-level two-mode generalization of the Jaynes- 
Cummings model. 



1. Introduction 

The definition of the canonical coherent states, CS for short, of the harmonic 
oscillator has been generalized in a number of ways and successfully applied in 
quantum physics for a long time. For various generalizations of this concept and 
their properties one may consult [SJ E El El . Hilbert spaces are the under- 
lying mathematical structure of several physical phenomena. In general CS form 
an overcomplete family of vectors in the separable Hilbert space of the physical 
problem. 

Let .ft be a separable Hilbert space with an orthonormal basis {0m}m=o ^^'^ *C 
be the complex plane. For z E T), an open subset of C, the states 



(1.1) |z)=AA(|z|)-^g 



are said to form a set of CS if 

(a) The states | z) are normalized, 

(b) The states | z) satisfy a resolution of the identity, that is 
(1.2) / \z){z\dfi^I 

where A/'d^l) is the normalization factor, {p(TO)}m=o a sequence of nonzero posi- 
tive real numbers, dfx is an appropriately chosen measure on J) and / is the identity 
operator on Sj. 

In a recent article jS], the labelling parameter z of (|1.1|) was replaced by an 
nx n matrix valued function of the form Z = A(r)e*^®*-'^-' , where A(r) and 0(fc) are 
n X n matrices and thereby a class of vector coherent states (VCS for short) was 
generated as n component vectors in a Hilbert space Sj = C" ® fti, where Sji is an 
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abstract separable Hilbert space. In order to get the normalization and resolution 
of the identity following conditions were imposed on the matrices. 

(1.3) e(fc) = e(fc)t, [^(r),e(fc)]=0, [A{r),A{r)^]^0, 

where f stands for the conjugate transpose of a matrix. The square bracket is equal 
to zero if the matrices commute. A further generalization of the same concept was 
given in |7], where the VCS were realized as infinite component vectors in a suitable 
Hilbert space and as special cases VCS were also given as finite component vectors. 
It is shown in [H] and |Zj that the VCS obtained with the complex representation 
of quaternions are very effective in describing two level atomic systems placed in a 
single mode electromagnetic field. In particular [7] argues that it is a good fit for 
the Jaynes-Cummings model. 

Gazeau and Klauder |3j proposed a class of temporally stable CS for Hamiltonians 
with non-degenerate discrete spectrum by introducing a generalization of the canon- 
ical coherent states. At present this class of CS is known as the Gazeau-Klauder 
CS. Following their method CS were derived for several quantum systems, see for 
example [S] |S] . As a further generalization of this concept multidimensional CS 
were introduced in [5] and as a special case of this generalization CS were obtained 
as a tensor product of two canonical CS in the following form 

(1.4) |.i,z,)=e-(l-l^+l-l^)/^^ Yl ^^l-i,«2). 

ni=0n2=0 V J- V ^ 

The result of can be considered as an extension to the results of 0] in a sense 
that it allows one to construct CS for a system with several degrees of freedom. For 
example, for a system with two degrees of freedom it can be written as follows: 

jni/2 

(1.5) I Ji, J2,7i,72) =M(Ji,J2)'^ V J--- e-'^^^"iAA2(J2,ni)-^ 

m VPi("i) 

VP2(ni,n2) 

where the labelling parameters and the Hilbert space were chosen appropriately for 
the system and p{ni), p{ni,n2), En^^, are associated with the spectrum E„j^„2 
of the system. The dependence of the first sum on the other severely restricts one 
to consider the states of Eq. (|1.5() as a tensor product of two CS. The detailed 
explanation can be found in 5 . This procedure is also used in |10) to obtain CS 
for a free magnetic Schrodinger operator. It should also be mentioned that in jll| 
Klauder has obtained CS for the hydrogen atom using a similar method of 
Motivated by these recent developments we generalize the results of [n|,[7] and |S] 
by introducing multi-matrix VCS in the form 



(1.6) I Ai,A2,...A,,j) ^M-i J2 E 



00 



. A "1 

\ ^ ^1 



ni— On2— Ut—O 

A "2 AKt 
-^2 ,J 



on a Hilbert space Sj = C" (8) [^l-^i^k], where A^'s are n x n matrix valued 
functions on appropriate domains, {0nj }5^^=o orthonormal basis of the sep- 

arable Hilbert space Sjj, j = l,2,...r and {pj(f^j)}i^=0' J — l:2,...r are positive 
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sequences of real numbers. A similar result with complex numbers zi, Z2, ■■, Zr may 
not be hard to obtain because the complex numbers commute with each other. For 
matrices, the non-commutativity restricts the construction severely. Furthermore, 
notice that 

I Ai{n),A2{r2),...Arirr),j) ^1 | ^2(^2), j) | Ar{rr),j). 

Next we discuss the physical motivation of the generalization (|1.6(l . It is well known 
that the most elementary model describing the interaction of two-level atoms with 
a single quantized field is the Jaynes-Cummings model. This model is exactly 
solvable in the rotating wave approximation, see |12[ I13L and the references 
listed therein. Suppose we have a diagonalizable Hamiltonian H for a two-level atom 
in a single mode cavity field with non-degenerate energies and wavefunctions 
ip'^, k = 1,2, m = 0, 1, ...,00. Let = E^ - E^, k = 1,2, m = 0, 1, ...c». 
Let R{m) — [diag(x^!, 2 and Z = diag(zi,Z2) be diagonal matrices, where 

x^! = x\...x'^ is the generalized factorial. Assume that the following vectors are 
normalized and satisfy a resolution of the identity, 

00 

(1.7) \Z,k)= U{Z)-i J2 R{m)Z"'xk 0m G C fc = 1, 2, 

rn—O 

where the wavefunctions V'm: fc = 1,2 are identified to the basis of ^ as 
4'm ■— Xk ^ 4'rm k — 1,2, 171 — 0,1, 2... The collection of vectors 

2 

|Z)-^Cfe|Z,fc) with |ci|2 + |c2|2 = 1 
k=l 

forms a set of CS for the diagonalized Hamiltonian Ho. Such an argument for a 
two- level system leading to the quaternionic VCS of ,£4 was given in [Jj. 

If one neglects losses, multi-mode multi-level generalizations of the Jaynes- Cum- 
mings model can be solved exactly using the exact solvability of the Jaynes- Cum- 
mings model |13l I15| . By assuming the existence of a solvable multi-mode multi- 
level system, we also construct VCS as CS of a multi-mode multi-level quantum 
model as an application to the multi-matrix VCS. In fact, an extension of p. 7(1 
to multi- level multi-mode systems requires multi- matrix VCS. For example, let us 
consider the following model. In the rotating wave approximation the Hamiltonian 
of a two-level atom in the two-mode field is given by jl5| 
(1.8) 

2 



i=i 



iJ = ^ f altti + ^ ) + x^o^a + .9 ( \/ a|aia^a2a|a^<T_ + a+aia2\J a\aiala2 



where H — 1, ai,al {i — 1,2) are the annihilation and creation operators of the 
iih mode of the field, uji (i — 1,2), cjq = ^1 + ^2 (i-e, the two-photon resonance 
process) are the frequencies of the zth mode of the radiation field and the atomic 
transition, 173, (t_ and (T+ are the pseudospin operators of the atom and g is the 
coupling constant. By diagonalizing H one can obtain the spectrum of H (see |15) 
and the references listed there): 

(1.9) E^:"^ =uoi{n + l)+ uj2{m + 1)± g(n + 1)(to + 1). 

For the zero coupling {g — 0) it can be easily seen that if ^ is an irrational 
number then the spectrum is non-degenerate. By appropriately identifying the wave 
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functions of H1.8|l as the basis of separable abstract Hilbert spaces, as an example 
of the general model, when g = we shall introduce CS for the Hamiltonian p.8|l 
in the following form: 
(1.10) 

OO oo 

I z, 3, k) - N{z, 3)-^ J2 Ri^y'^z^ J2 "*)"^3'"xfc (8) fc = 1, 2, 

n— m— 

where Z, 3, Rin), R(ri, m) are 2x2 diagonal matrices. In addition to this example, 
we shall also present VCS for tensored Jaynes-Cummings type Hamiltonians. 
In order to avoid technicalities, we carry out our general construction with two 
n X n matrices A and B. Later using the real matrix representation of quaternions 
and octonions we demonstrate it for several matrices. We also generate VCS in the 
above form with matrix valued /9j(nj)'s. Further, we discuss VCS of type (|1.5|) . We 
shall briefly discuss generalized oscillator algebras associated to the VCS. Finally, 
we introduce VCS for multi-level multi-mode quantum systems. 

2. General construction 

Let be two separable Hilbert spaces and C be the complex plane. Let 

{4'm}m=0j {V'ij/^o ^^'^ {X"'}5'=o be orthonormal bases of the Hilbert spaces in the 
respective order and of C". Denote S) ~ C" (g) i^i (X) io2, where (g) stands for the 
tensor product. Then the set of vectors {x-' 'S) 4>m 'S) V'/}i,m,; serve as an orthonormal 
basis of Sj. 

Let K, K' be measure spaces with probability measures dK and dK\ respectively, 
and i?, S" be a second pair of measure spaces with measures dR and dS, respectively. 
For {r,k,C) e Rx K X [0, 27r) and (s, k', rf) e S x K' x [0, 27r), let 

(2.1) A = A(r)e*'^®W, and ^ = B{s)e''^^^'''\ 

where A(r), 0(fc), _B(s), A(fc') are measurable n x n matrix valued functions with 
the following properties (assumed to hold almost everywhere with respect to the 
corresponding measure(s)): 

(2.2) Q{k)^Q{k)\ [A(r),A(r)t]=0, [A{r),e{k)] = 

(2.3) A{k')^A{k')\ [B{s),B{sy=Q, [B{s), A{k')] = Q 

(2.4) [B{s),A{r)]^Q, [A{k'), A{r)] ^ 0, [S(s), A(r)t] = 0, [e(fc), S(s)] = 0. 
Let us recall that 

(2.5) [Bis), Air)] = 0, [A(fc'), ^(01 = 0, [Bis), A(r)t] = 0, [e(fc), Bis)] = 
=^ [Bis)\Air)^]^0, [A(fc'),A(r)t]=0, [S(s)t, A(r)] = 0, [e(fc), B(s)t] = 
Let V — RxSxKxK'x [0, 27r)x [0, 27r) and define the measure rf^(r, s, k, k' , (, rf) = 

dRdSdKdK' dC^dr] on it. For each pair of matrices A, B we define multi-matrix vec- 
tor coherent states (MVCS) as follows. 

oo oo 3^ 

(2.6) I A,B,j)=AA(r,s)-5 ^^^_^==X^"®0™®V'i, 

where the normalization factor A/'(r, s) and the positive sequences pim),\il) have 
to be chosen so that the MVCS are normalized and give a resolution of identity. 
In the sequel we use the following notation: For a,n n x n matrix A, we denote 
|A| = [AA^h. 
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In order to normalize the states H2.6|l and to obtain a resolution of the identity 
we impose the following assumptions on the matrices A and B: Let A and B be 
matrices of the form (|2.1|l satisfying conditions (|2.2|) . H2.3f) and l|2.4|l . Suppose that 
Q(k) and A(fc') have non-zero integer spectrum. With these conditions the vectors 
(12. 6() are normalized to unity as 

n 

(2.7) ^(A,B,j|A,B,j)=l 

i=i 

with the normalization factor 

m=0 1=0 ' ^ ' 

In order to see this let us start from 12.7|l and establish 12. 8|) in the following way: 

n n oo oo oo oo 

^(A, B, J I A, B, = U{r, s)-'Y.T.T.T.T. ^ ^( Ur/Ur ^ 

j = l j = l m=0 /=0 .=0 r=0 VpMp(^^)^(OA(t) 

n OO OO OO OO ^ 



t m=0 i=0 .=Q r=Q VPMpWMOAM 

X(A™B'X^' I A'^B-x')cA^ra I ^.)S,A^I I V'r)fl. 



n OO OO 



because 



(A^B'x-'' I A^B'x^k" = (B'tA"tA"B'x-'' | x')c'^ 
(B'tA(r)"tA(r)"B'x-'' | x'>C" [by (E3] 
(B(s)'U(r)™tA(r)™S(s)'x^' | x'')c" [by lOl and (lOl l 



Therefore, we have (|2.8|l . For the resolution of the identity, we have 
V / I A,B,j)(A,B,j I d/^(r,s,/c,fc',C,??) 

noooooooo ^ p p p p p27r i^lix 



1 



§ £ S S S Vp("^)p(^)A(0 A(r) }b. Js Jk Jk' Jo Jo AA(r,s) 
<8) I (/>m) I «> I -0;) (V't I s, A:, fc', C, »7) 
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DC OO OO OO ^ 

^^^^ — , , , , , , 

^0 ^0 ^0 Vpim)piiy)Xm{T) Jr Js Jk Jk' Jo Jo AA(r,s) 

xA(r)"'e""'^®('=)B(s)'e*''"^('''')B(s)^tg-^Tr,A(fc')^(^)i^te-''Ke(fc) 
«> I (f>m) I «> I -0;) (V'r I dfi{r, s, k, k' , (, rf) 

n 

because | x'){x^ H nx n identity matrix 
i=i 

47r2 f f f f 1 

xA(r)™B(s)'B(s)'tA(r)™t,^ | 0™)(</)™ | ® | I dRdSdKdK' 



by l|221), (ESll, JZl and ^ e'('"-'')^®('=)dC = | 



if m ^ v 
27r if TO = J/ 



CXD OO 



[because dK, dK' are probability measures, and by H2.4|l ] 

provided that the measures dR and dS* satisfy 

(2.9) 47r2 / / H!ffl:?ifl^di?d5 = p(TO)A(Z)I„. 

iflis A/(r,s) 

Thereby, the situation can be summarized as follows: The set of states (|2.6f) forms 
a set of MVCS if N{r, s) of (|2.8|l is finite and the condition (|2.9|) is satisfied. In this 
case the states are normalized as ()2.7f) and a resolution of the identity is obtained 

as 

n „ 

V / I A,B,j)(A,B,j I d/^(r,s,/c,fc',C,?7) =I„0/f„ 

In the following let us see some illustrative examples of the main construction. 

2.1. Multi-quaternionic VCS with complex representation. As an example 
of our general construction, in this section we derive MVCS using the complex 
representation of quaternions by 2 x 2 matrices. Using the basis matrices, 

10\. f Q i \ . f Q ~l\ . f i 



^0 = 1 1 ; ' = 1, . ; ' ^[i J' '"'^yo 

where ai,a2 and (T3 are the usual Pauli matrices, a general quaternion is written 
as 

q — XQffo + ix- a 
with aiQ G M, x = (xi, 0:2, 2:3) € R'^ and ct = (ai, — (T2, cs). Thus, 

xq + ix3 —X2 + ixi 

X2 + ixi Xq - iX3 

It is convenient to introduce the polar coordinates: 
a;o = rcos6', = r sin6'sin(/)C0S V', a;2 = r- sin 6* sine/) sin -0, X3 = r sin 61 cos , 



(2.10) 
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where r e [0, oo), <f) £ [0, tt] and Ojifj e [0, 2tt). In terms of these, 

(2.11) q = A(r)e*^'"(") 

where 



/^-,r,x ^ / \ / COS0 sin0e"^ \ , ,^,9 

(2.12) A(r) ^ rao , ^(n) = _ ^^^^ j and a(n)^ = a, . 

We denote the field of quaternions by H. Now let us take two quaternions qi, q2 as 
follows. 

(2.13) qi = A(r)e^^i"("i\ qa - B(s)e*^^"("^\ 
where 

., , ,^ \ ( cos(?!)i siniTiie*'''^ \ „, ^ 

Ayr) ^ rao, a[ni) ^ [ . , , «(s) = scfq, 



and CT(ri2) = 



cos 02 sin(/)2e''^^ 
sin (/)2 e ~ ''^^ — cos 02 



The matrices A{r), B{s), cr(rii) and (7(^2) satisfy the conditions H2.2|l - 12.4|) and the 
eigenvalues of cr(rij); j = 1,2 are —1 and 1. Let {^m}m=0'{^ iJl^o ^'^ orthonormal 
bases of abstract Hilbert spaces i5i,^2 respectively, and x^,x^ be an orthonormal 
basis of C^. Let V ~ [0, 00) x [0, tt) x [0, 27r) x [0, 27r). Let us take the measure to 
be 

W(r, s) 

d/x (r, s , 01 , ■01 , ^1 , 02 , "02 , ^'2 ) = \ rdrsds (sin 0i )c?0i dipi dOi (sin 02 ) d02 d0'2 ^6*2 

Idtt-^ 

on V X V, where W{r, s) is a density function to be determined. With the above 
setup we have a set of MVCS, 

00 CXD ^ ; 

(2.14) |qi,q2,j)=AA(r,s)-^^5]-E^-%=X^'®'i>™®vl/,, j ^ 1,2. 

From (|2.12l) . 12.13|l and (|2.8|l the normalization factor can be obtained as 

00 00 

(2.15) AA(r,.) = 2^^ 
Since 



p{m)\{iy 
m=o ;=o '^^ ^ w 



''\'('"-'')^^-("^)sin0id0idM^i = / ^''^^ if m = i. 



"'0 



if TO 7^ 



and a similar integral holds when the index 1 changes to 2, from (|2.12|l . (|2.13l) and 
(|2.9|) . following the argument of Section [3 when W{r,s) satisfies 

j-oa /-oo TIT/ N 2m+l 2/+I 

(2.16) 47r2 / / "^^"^'^r. / drds = p{m)\{l). 

Ja Jq -^{r, s) 

we can readily obtain the resolution of the identity 



(2.17) / I qi,q2,j> M/(r,s) (qi,q2,j M/^ = 12® 

JVxV 

In the particular case p{m) = T{m + 1) and X{1) = T{1 + 1) we have 

AA(r,s) = 2e'''e"' = 2e'^'+^' 
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and W{r, s) — ^ satisfies (|2.16f) . which can be seen with the integral representation 
of the gamma function. 

Remark 2.1. • One could solve the moment problem (|2.1t)|) for various choices of 
p(m) and X{1) using the methods demonstrated in [16j and thereby create several 
classes of multi-quaternionic MVCS. 

• While keeping conditions (|2.2|) . 12.3|) . if we replace conditions (|2.4|) by 

(2.18) ^(r)A(r)t A(r)U(r) = /(r)I„ and B{s)B{s)^ ^ B{s)^B{s) ^ 5(s)I„, 

where /, g are functions of r and s, we can carry out the construction in a somewhat 
similar way. We demonstrate it through the following section. 

2.2. An extension of the quaternion case. In the case of quaternions we had 
a more convenient form for the matrix A{r) in which the matrix A(r) is simply 
r times the identity. Here we will establish MVCS with a little more complicated 
j4(r). This example is, in fact, a generalization of the quaternion case. Let 

" rh -sh ■ 



6 = 6(ni,n2, 



cr(ni)sin6' ia{n2) cos9 
—iaiji,) cos 9 a{ni)si'a9 



where ni — (nii,ni2,ni3) and ~ («-2i, ^^22, '^■23) are unit perpendicular vectors 
and a{nj) = (riji, n-,-2, njs) ■ (cti, (T2, (72). Let 

Bi = v4(ri,si)e<i«i and B2 = A{r2, S2)e"^'^\ 

where Qi = Oi(n']^\n'^\9i), 82 = 02(jif\n^\92)- Through straightforward 
calculations we can see that A{ri , si), A{r2, S2), Oi and 02 satisfy conditions (|2.2|) . 12.3|l 
and (|TTH|) . Letri,si,r2,S2 € [0, 00), Ci, 6*1, C2, 6*2 e [0, 27r) and n^^\ \ n^^' e 

il, where fl is an open subset of M.^. For the measure let us take — dvdndVl^ 
where dv = riS\r2S2dr\dsidr2ds2, dn — -^dC,idC,2 and dVl is the probability mea- 
sure on X X [0, 27r) x [0, 27r). Let 

V = x [0, 27r) X [0, 27r) X X r2 X [0, 2tt) x [0, 27r). 

With the above choices and with the usual notations wc form a set of MVCS: 
(2.19) 

\Bi,B2,]) =AA(ri,si,r2,S2)-^ V V ^^^=|=X^®0m'55V'i, J = 1,2,3,4. 

m=o 1=0 VPi'^) V^(0 

Since 

BfBT = [rl + s?)"l4, and S^tg?, = (rj + sl)%, 

and 

{BTB'2X' I BTB'2X') = {B'2\BfBTB^2X' I x') = ^{rl + s?)™(r: 



we have 
(2.20) 



4 °° °o /„2 I „2\m („2 , „2y 

5:(S„S2,, I B„B2,j)=mn,s„r2,S2r EE ,1^) Am 
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Thereby the normahzation factor is given by 



m=0 1=0 ' W 



Further, since 

4 



we have 



4 /■ 

V / W{r^,s^,r2,S2)\Bl,B2,j){Bl,B2,]\d^l 
i=i 



m=Ol=0 P^'^'>^^^^ ■^^ -^0 -^0 -^0 A^(?-l,Sl,7-2,S2) 

«) I (/'m)((/'m I «) I ll>l){ll>l I . 

Thus a resolution of the identity is satisfied provided that there exists a density 
W{ri,si,r2,S2) satisfying 

^oo /'oo /'oo roo TJ// \/ 2 i 2\m t 2 i 2\/ 



Ar/ \ '<^^ = p{m)X{l). 

JO JO JO N{ri,si,r2,S2) 

For the particular case, p{m) = m! and — U the normalization factor can be 
easily deduced from H2.2()|l . 

AA(ri,si,r2,S2) = 4e(''i+^?+'-'+^') 

and a resolution of the identity is obtained with 

W{ri,si,r2,S2) = W{ri,si)W{r2,S2) 

where 

W{ri,si)^— — 2> W^(''2,S2) = — 2' 
rf + si ri + s\ 

Which can be seen in the following way: Since in general we have 



M{B) 



-srdsdr 



Jo 

roo roc ^g(^2^^2)m 

-srdsdr 



4(r2 + s2)e'-^+s^ 

oo r°o 4(^2 ^ ^2)m-l 

^ — ^ srdsdr 



10 Jo 

m — 1 



i— ^ / 1 \ POQ poo 



m — 1 



3=0 
m — 1 



3=0 



J 



E ( "'j ^ )r(m-j)r(j + l)=m!. 



OO / /"OO 



Vo 
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we obtain 



4 

E 



/ Wiri,Si,r2,S2) I Bi,B2,j){Bi,B2,j \ = I4 » ® . 
Jv 



Remark 2.2. For r, s G [0, 00) and 6', 77 G [0, 2n), if we take our matrices A and B 
as, 

(2.21) A = e*^A(r) and B = e*''B(s) 
with conditions, 

(2.22) A(r)A(r)t = yl(r)U(r) = /(r)I„ and B(s)B(.s)t = B{s)'<B{s) = g(s)I„, 

we can carry out the above construction with conditions (|2.2f) . H2.3|) and (|2.4I) re- 
placed by the single condition H2.22|l . The condition 12.22|l is satisfied by any real 
Clifford type matrix |17l [7] . Following the construction given in Section |21 one can 
construct MVCS with any n x n Clifford matrix. However, since an explicit real 
matrix representation is known only for quaternions and octonions '17' , in the fol- 
lowing sections we demonstrate it with the real matrix representation of quaternions 
and octonions. 

2.3. Multi-quaternionic VCS with real representation. Here we present quater- 
nionic MVCS with the real matrix representation of quaternions without imposing 
any conditions on the matrices. Let 

H = {q' = ao + flii + a2j + a^k \ — = fc^ = —1, ijk = —1, ao, ai, a2, as G M} 

be the real quaternion division algebra. It is known that H is algebraically isomor- 
phic to the real matrix algebra 



M={q' = 



I ao -oi 

ai ao 

0-2 as 

V «3 -02 



-a2 
-as 

ao 

ai 



-as ^ 

a2 
-ai 

ao / 



ao, ai, a2, as G 



For a detailed explanation see ^7] and the references therein. For z = e'^ G and 
q' G M, let q(z) — zq' , then we have 

q(z)q(z)t = q(z)tq(z) = (a^ + aj + 4 + al)h = \q\% = \q'\% 

where |q'| is the norm of the quaternion q' and I4 is the 4x4 identity matrix. 
Let q'(a) denotes the matrix representation of q' with matrix elements a^s. Let us 
consider the Hilbert space (8) i3i f)2- The quaternion norm can be considered 
as a continuous function. 



qt-i- 1. 



Thus, we make the following identification: |qi| — t, |q2| — s. Since zi = e ^ and 
Z2 = e**2, set the measure dfj,{t, s, 6*1, ^2) = d9id92tdtsds owV — [0, 27r) x [0, 27r) x 
R+ X K+. With the notations of the previous section let us define a set of multi- 
matrix VCS as follows: 

For j = 1,2,3,4 the following set of vectors forms a set of MVCS. 



00 oc 
1 \ — > \ — > 



(2.23) I qi(zi,a),q2(z2,6),j) ^ M {\q,l\q2\)-- ^ 



qi" 



1=0 vpM VW)' 
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Let us find the normalization factor. Since 
4 



^{qi{zi,a),q2{z2,b),j \ qi(zi, a), q2(z2, j) 

1 



4 oc oo 

A^dqil, Iq2|)-^ E E E 7;7;;;wn > i^^^x^) 



j=l m=0 1=0 

|qir"'iq2 



m=0 i=0 '^^ J \ / 

the normahzation factor is given by 

°° °° |„ |2m|„ |2Z 

,2.24, A.,|,.|.te|,=4EE^blr- 



Next we obtain a resolution of the identity. Since 

E / I qi(^i'«)'q2(^2,6), j>(qi(2;i,a),q2(2;2,fe),j MM(i,s,^'i,6'2) 

4^2 .oo .oo^(^^^2™+1^2/+l ,,,,,, 

^o^pM^lOio Jo JV{t,s) 
provided that there exits a density function W{t, s) such that 

poo poo Txr/, \j.2m+l 2/+1 

(2.25) An' / / dtd.^p(m)A(0. 

Jo Jo A/ (t, s) 

Remark 2.3. • For the choice p{m) = to! and A(/) = /! from H2.24|l we have 
7V(|qi|, Iqsl) = 46*'+^*' and W{s,t) = satisfies 

• One could easily notice that the whole procedure depends on the following: 

(2.26) (qrq^X^' I c^Tq^x') = ((qrq2)^qrq2X^' I x') = IqiP^lqsP' and 



4 

(2.27) E N^l2X^Kl^l2X^' I 



The identities H2.26|l and 12.27|l holds for any number of quaternions, that is, qiq2 
can be replaced by qi....q7., therefore the procedure can be extended to have MVCS 
on the Hilbert space C" (8) [^r=i-^r], where S)r, r = 1, ...,t are separable Hilbert 
spaces. 

2.4. Multi-Octonionic VCS. Let O denote the octonion algebra over the real 
number field R. In |17| it was shown that any a € O has a left matrix representation 
a; (a) and a right matrix representation iy{a) and these representations were given 



qi q2 




x'){x' 



(qrq2)^ 



= qi q2q2 qi = lqil Iq2| u- 
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respectively by 



uj{a) = 



v{a) = 



/ 

/ Oo 


-ai 


-a2 


-as 


— a4 


-as 


-ae 


-ay 


\ 


ai 


ao 


-as 


a2 


-as 


a4 


ay 


-ae 




a2 


as 


ao 


-ai 


-ae 


-ay 


04 


as 




as 


-a2 


ai 


ao 


-ay 


ae 


-as 


a4 




04 


as 


ae 


ay 


ao 


-ai 


-a2 


-as 




as 


— a4 


a? 


-ae 


ai 


ao 


as 


-a2 




ae 


-ay 


— a4 


as 


a2 


-as 


ao 


ai 




V ar 


ae 


-as 


—04 


as 


a2 


-ai 


ao 


/ 


( ao 


-ai 


-a2 


-as 


—04 


—as 


-ae 


-ay 


\ 


ai 


ao 


as 


-02 


as 


— a4 


-ay 


ae 




02 


-as 


ao 


ai 


ae 


ay 


—04 


-as 




03 


a2 


-ai 


ao 


ay 


-ae 


as 


— a4 




04 


—as 


-ae 


-ay 


ao 


ai 


a2 


as 




as 


04 


-ay 


ae 


-fli 


ao 


-as 


a2 




ae 


a? 


04 


-as 


-02 


as 


ao 


-ai 




\ a? 


-ae 


as 


a4 


-as 


-a2 


ai 


ao 


/ 



Let z = e*^ e S'l and 

aj(a, z) = ZLo{a) and :/(a, z) = zv{a). 

Then 

ti;(a, z)a;(a, z)^ = Lj{a^ z)'^ Ljj{a, z) — i'{a, z)v{a, z)^ = via, z)'^ v{a, z) 

= + a^ + a2 + ag + a| + ag + Og + a^ = ||a|pl8. 

Thus, identities similar to H2.26|l and (|2.27|l can be obtained for a;(a, z) and z^(a, z). 
Now, with previous notations we can have MVCS as 

I uj{ai, zi),a;(a2,Z2), t^(ar, z^-), j), 
I z^(ai,zi),i/(a2,Z2),...,i^(ar,z^),j), or 
I z^(ai, zi), a;(a2, Z2), i/(as, zs), 0^(04, Z4)..., ^{ar, Zr)J)- 
In the last set of states we can mix w(a, z) and i>{a, z) in any order. 

3. Multi-matrix VCS with matrix p{m) 

Here we consider a class of MVCS in the following form 
(3.1) 

00 00 

\ Zi,Z2,..,Zr,j) =M-^ - J2 Rl{ni)Z^\..Rrinr)Z!^-x' ^ (l^n, ^ - ^n^, 

711—O n^— 

where Rk{nk) and Z^ are nx n matrices for all fc = 1, r. In order to get a 
normalization and a resolution of identity, we need to take the ZkS in a certain way 
and we have to impose conditions on all the matrices. For example, for /c = 1, ...,r 
let Ck be a real n x n Clifford matrix, Zk = ZkCk with Zk G and Rk{nk) be a 
n X n matrix satisfying 

(3.2) Rk{nk)Rkink)^ = RkiukY Rk{nk) = fk{nk)\, 

where /fc (n^ ) is a real valued function of rifc . Then following the argument of Section 
El one can construct MVCS. For simplicity, we work it out in detail with the real 
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matrix representation of quaternions. Let us take qfc to be as in Section 12.31 for all 
fc = 1, 2, .., r and 

(3.3) Zfe = qfe(^?fe) = qfee'^'=. 

Further, let us take Rk{nk) such that 



(3.4) Rkink)Rkink)^ =Rkink)^Rk{nk) = fk{nk)h = ^Ii, A: = 1, 2, r. 



For example one such possible Rk{nk) is 

1 / cos(a;fc)l2 — sin(a;/c)l2 



Rk{nk) 



Ukl V sin(xfc)l2 cos(xfc)I 



2 



where xi, Xr are fixed. 

For fc = 1, T let Zk be as in (|3.3|) and Rk{nk) be as in (|3.4|l then the vectors (|3.1|) 
form a set of MVCS. Since 



{R,inMeir\..Rr{nr)qAOrr-x' I i?l(nl)ql(0l)"^..i?r(nr)qr(^?r)"-x''> 
= (q,(0O"^^^r(rv)t...qi(^?i)"^ti^^(„j)t^^(„^)q^(0^)ni...^^(„^)q^(0^)«.^j- | ^J) 

|qi|2"i...|q,|2". . 



ni!...nT-! 



-(x^ I X' 



and 



I i?i(ni)qi(0i)"^..i?.(n.)qr(0r)"-X^)(i?i(Zi)qi(^i)'^..i?r(Zr)qr(0r)'^x'' 



xql-^Rr{lry...q['^Ri{h)^ 
e^("i-'i)«i...e^("--'-)«-i?i(ni)q!^i...i?^(n^)q;!- x q^-t^^(^^)t_ .q'^it^^(^^)t 



we have the normalization factor 



4 °o °° \„ |2ni |„ |2n^ 

(3.5) M = J:J:-E ''^^L I X^) - 4el-l^--l-l^ 
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Set dfi — j:iy\qi\...\qr\d\qi\...d\qr\d9i...d9r- With this measure a resolution of the 
identity is obtained as follows: 



E / - / / - / |qi(^i),-,q.(^r),j>(qi(0i),...,qr(er)j |dM 
^^-^ Jo Jo Jo Jo 

OO OO OO OO 

E-EE-E 



OO OO OO cxj „oo /.OO p2tt p2tt 



1 

ni=0 n,=0/i=0 1^=0-"' " " " " " " 

® I I I <j)r){(t)r I 



= E-E2^/ ■■■/ e-l.^|-...e-l..|-|q,|...|q,Hq,|...rf|q,| 

(g) I I I <j)r){(t)T I 

Remark 3.1. Usually the p(to) of (|1.1|) is a positive sequence of real numbers and 
in getting a resolution of the identity we end up at solving a moment problem with 
positive moments. When we replace the p{m) by a matrix, one could expect to 
replace it by a positive definite matrix. The matrices we have used are not positive 
definite. But one can observe that, at the end of the calculations, we end up with 
a moment problem with positive moments. In this sense, it really does not matter 
with which matrix we start with but what does matter is the final moment problem. 

4. Summations depending on one another 

So far we have demonstrated several classes of MVCS as a generalization of H1.4|l . 
In this section we present MVCS as a generalization to the states H1.5|l . where the 
summations depend one another. First we give a general construction and then we 
discuss an example using the complex representation of quaternions. 
Let A, B be matrices as in Section |21 with all the assumptions of Section[21on them. 
Then the following set of vectors forms a set of MVCS. 

(4.1) I A, B,j) =AAi (r,s)-^ V ^^^ J^2 (g , m) ' ^ V —==x' ^c^m^i'i, 

where /, g are some functions of m and I respectively. Following the steps of Section 
121 we can get the normalization factor 

(4.2) M(r,.) = 5]^ ^'1^^"^ ""^'^ I 



Again by following the steps of Section |21 we can have the resolution of the identity 
V / \A,B,j){A,B,j\dp{r,s,k,k\C,v) 
^ ^ 47^2 r f |A|2/(™)|BP9« 

E E 7r7;;;w;;rn ll .)Limf ' ' ^ ' ' 

xAi(r)A2(r, s)dRdS 



„=o /=o PiMP2(m,0 JjiJs Afi{r,s)Af2{s,m) 
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provided that there are densities Ai(r) and A2(s,to) such that 

(4.3) 47r2 / / ' . -\iir)\2{s,m)dRdS ^ piim)p2im,l)K. 

JrJs s)7V2(s,r7i) 

where R, S and dR, dS are as in Section [5] 

Now let us see an iUustrative example with the complex representation of quater- 
nions. Take the quaternions as in H2.11|l . Let 

01 = {q = yl(r)e*'^'^(^i) : rG[0,oo)}, 02 = {p = A{s)e''"'^"'^ : se[0,l)}. 

The variable s is restricted on the set [0,1), the rest of the variables obey the 
conditions of SectionEH For q G (5i, p G ©2, /(m) = ^,g{l) = 5 with 

f m + I \ ^ 

A = q, B = p, p2{m,l)=l ^ I , /9i(m) = m! and j = 1, 2 

let us consider the states H4.1II . From (|4.2|l and the properties of quaternions we 
obtain 



00 oc 



(4.4, = -^'^^MiJli^ 

m=0 1=0 ^ ^ 



where 

oc 

(4.5) N2{s,m) = Y,[ )5" = (l-5)-("+i) 



1=0 

Therefore we get 



(4.6) 



M,{r, s)=2Y: ^ > , = (1 - ^)e'^(^-^) 



m=0 



For a resolution of the identity (|4.3|) reduces to 

(4.7) 4^2/ / (r) A2 (r, s, m)drds = m! ( ^ 



Since 

(4.8) m / s'(l-s)™-ids = 







m + Z 
I 



we take 



A2(r,s,m) = — — -e 

27r(l — s) 



Thereby from (|4.7|) we only have to solve 

/"OO 

(4.9) 2n r™e"''Ai(r)dr = m!. 



By the definition of the Gamma function, we obtain (|4.9|l when Ai(r) = Thus 
we have a resolution of the identity. 
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5. Generalized oscillator algebra 

For the states , a general way of defining an associated oscillator algebra is 
as follows pP: Let 

- - ^^""^ Vrn > 1, and xqI = 1 



p{m — 1) 

then p{m) = Xm- A set of operators, namely annihilation, creation and number 
operators on the basis vectors {4>m} is defined as 

a0m=\/^0m-l, = V^m+10m+l, Ntpm^Xm^'m- 

The states satisfy the relation a\ z) = z \ z), that is, the CS are eigenvectors 
of the operator a. These three operators together with the identity operator, under 
the commutator bracket, generate a Lie algebra, which is the so-called generalized 
oscillator algebra. In analogy with the above case, for the states (|1.6f) let us take 

- - ^^"^"^^ .o! = l; V, = l,...,r 



Let us define a similar set of operators for the basis vectors {x^ (pm (8) ■■• 
as follows: 

(5.1) Ax' «) (j)ni ... «) (pn^ = y/XniXn2---Xn-,X^ ^ Ki-l ® ■■■ ® 4>n^-l 

A^X' ® 4>ni ® ••• = ^/X„^ + lXn^^+i...Xn^ + ix' (g) (f>m+l ■•• » 0n,+l 

Nx' ® (j)ni •■■ = XniXn2---Xn^x' ^ ■■• (t>n^- 

If 

(5.2) Xrij — Xrij-i = c, a constant = 1, ...,t, VjIj 
we get 

[A, A^ = cl, [N, A] = -cA, [N, A^] = cA^ 

and the algebra is isomorphic to the Wycl-Hcisenberg algebra. Since Aj 's do not 
commute, in general, we have 

A I Ai, As, .-ArJ) ^ AiA2...A^ | Ai, Aa, ...A^, j). 

Observe that if the matrices Ai,A2,...,At are replaced by complex variables 
zi,...,Zt and the basis is replaced by ^ ■■■ ^ <t>n^ then the states H1.6f) can 
be considered as multi-boson states with no interaction between bosons. In such a 
case, when (|5.2|) is satisfied, if we define a set of operators similar to the above case 
we can have 

(5.3) A I Zi,...,Zr,j) = Zi...Zr \ Zi , . . . , , j) . 

For the states (|1.6|l . if we assume that the matrices satisfy [Aj, A^] = for all pairs 
j ^ k we can have a relation similar to (|5.3I) . Such an assumption will be very 
strong. 

In order to have an annihilation operator so that the MVCS | Ai, A2, ...A,-, j) as an 
eigenvector of it, it has to be defined in such a way that the action of the operator 
A affects only the first component of the vector, that is, 

(5.4) Ax-' ® (pm (Xi ... <E) (pnr = Vx^X^ (pm-i (pn^ <^ ■■■ <^ (pn^- 
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In this case the other two operators take the form 

(5.5) ^t^-' (g) = ^Xn^ + ix' ® 0ni + l ® 0«2 ® ■■■ ® <t>n^ 

(5.6) Nx' ® 4>ni® 4>n^ = XniX-' </>ni (8)0„2 (8) ... (8) 0„^. 

Under the definition H5.4() we have 

(5.7) A I Ai,A2,...A^,j) Ai | Ai, A2, ...A^, j) 

and if we assume (|5.2|l we again get an algebra isomorphic to the Wyel-Heisenberg 
algebra. 

For the states of Section 0] one can define the operators by the same relations of 
(|5.4|) , H5.5|) , H5.6|l and get a relation similar to (|5.7|) . If we attempt to define operators 
as ()5.1|l we end at a situation similar to multi-boson case with interaction between 
bosons. 

6. A PHYSICAL MODEL 

In the following we adapt the above methods to construct a class of CS for a 
Hamiltonian, which is a tensor product of r multi-level quantum Hamiltonians. 
In the sequel we also derive CS for tensored Jaynes-Cummings models and for a 
two-mode two-level generalization of the Jaynes-Cummings model. 

For each j = 1, 2, r, let be a Hamiltonian for an n-level atomic system. Let 
-^k nij ^'^'^ ^kmj ,k — 1, ...,n\ TTij = 1, 2, oo bc the eigenenergies and eigensolu- 
tions of the Hamiltonian Hj for j = 1, ...,r. Suppose that there is no degeneracy 
and the energies are ordered as follows. 

< < < El^ < fc = 1, 2, n; j = 1, 2, r. 

For each j = 1, 2, r suppose also that the eigensolutions 4 , fc = 1, n; rrij = 

l,2,...,cx) form an orthonormal basis of the quantum Hilbert space ^gjv/ of the 
Hamiltonian Hj. Let us rearrange the energies as follows: 

4,m, ^H,m,^H,0' fc = l, TOj = 0,l,...,oo; j = l,2,...,r. 

Thereby we have 

= 4,0 < 4,1 < 4,2 < fc = 1: 2, n; j = 1, 2, r. 

Let be abstract separable Hilbert spaces with orthonormal bases {"^mj }mj=oi J — 

1, 2, r and {xk}k=i be the natural basis of C". Let 

r r r 

H^^Hj, S)qm^(^%m^ and i3-(g)%. 

Then _ff is a Hamiltonian describing n level tensored systems with the state Hilbert 
space SjQM- The eigenenergies and the eigensolutions of H can be written as 

r r 
<Pk,m^,...,m^ = (8)<^im. ^^'^ Ek,m^,...,m^ = W_K,-mp 

where fc = l,...,n. Observe that the energies Ek^mi,...,mr. of the Hamiltonian H 
can be degenerate. However, it will not affect the CS construction. Since we sum 
over all the indices of the energies and the energies of the individual Hamiltonians 
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Hj, j = l,...,r are non-degenerate, the degeneracy of Ek,mi,...,mr. will not cause 
any problem in obtaining a resolution of the identity. Define 

n oo oo 

V:S^QM >C"(g)i3 by ^ = XI ' ' ' I Xfe l/'mi,...,m,)(<?!'fc,mi,...,m J, 

fc=lmi=0 mr=0 

where 



^mi,...,m, = QPV'm^- 

Then V is a unitary operator with 

<^/c,mi,...,mr ^ Xk ^ V'mivjTOr 

for all fc = 1, n and mj = 0, 1, 2, j = 1, n. With the above setup we have 

VHV-^ =Hd= diag(ffP, Hi>, if^) 

a diagonal n x n matrix, where H^, k = 1,2, ...,n are self-adjoint operators on Sj 
with 

-^fe V'mi,...,m^ = -£'fe,mi,...,mrV'mi,...,mr5 = 1, 2, 77.; 771^=0,1,2,... 

Let S C and 0^ = r^e^*^, j = 1, 2, r; A; = 1, 2, n. Let 

= diag(2;^ , 0^, 4) and pj{mj) = diag(ei_„ . !, 4 !, 4,m,- ')> 
where sj.^^. ! = Sk,i^k,2 ' ' ' ^k,mj ■ es-ch j and k, suppose that the series 

™3=0 y^k,mj - 

is convergent and the radius of convergence is i^. Let 

D, = {{zi,zi,...,zi)GCx...xC I 141 <iU- 
Suppose that there are densities A^(4) such that 

(6.2) _^^'(4)2™Ai(4)d4=4„^.!; fc = l,2,...,n; j = l,2,...,r. 
Set 

^ r n 

d„{Z,, Z^) = —— n n Kiri)dridei 

^ ' j=i k=i 
on Di X ... X Dr- With this setup the set of vectors 
(6.3) 



00 00 



i=i 



\z,,Z2,...,Zr,k)=^^{Zu...,Zr)--^ X E ••• E 

mi— m2— m,,— 

■ Xk i^lni ® -^112 ® ■ ■ ■ ® 

forms a set of MVCS, when the states are normalized as 



^^(Zi, Z2, Zr, k I .^1, ^2, Zr, k) = 1, 



k=l 
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and the normahzation factor is given by 



|2mj 



mi— 0m2— m,,— Vi~^ 

With the argument of Section |2] and with H6.1() it can be seen that 

Af{Zi,...,Zr) < oo. 

Again with the argument of Section |21 and with (|6.2|l one can prove the resohition 
of the identity in a straightforward way as follows: 

n „ 

/ \ Zi,Z2,...,Zr,k){Zi,Z2,...,Zr,k\ 

(0-4) i^^^JDiXD2X...XD,. 

■jV{Zi,...,Zr)dn{Zi,...,Zr) =I„«)/^i Isj,.- 

Remark 6.1. Using the technique of jZj the diagonal nature of the VCS H6.3|l can 
be removed. For this let dv be the normalized invariant measure of SU (n) and let 
U e SU{n). In (|6.3|l let us replace 

n P,K)-^C^ by [/ n P,(m,r^Z7^C/t 
j=i i=i 

then a resolution of the identity can be obtained with the measure d^{Zi, Zn)di>{U) 
and the normalization factor remains the same. For a detail explanation we refer 

Remark 6.2. Let qi and q2 be two quaternions as in Section 12.31 Then it is 
straightforward to see that the eigenvalues of qi are zi = oq + ib and zi each 
with multiplicity 2, where b = y^j^ + a| + a|. Furthermore, there are orthonor- 
mal eigenvectors X^^^"* and Xj^^"*; J — l?^. Let Ui be unitary matrix such that 
qi = [/idiag(zi, zi, zi, ~ UiZiU\. Let the corresponding counterparts of q2 

be Z2 and U2- In 16.3|l one could take k — A, r = 2, Zi = diag(zi, zi, zi, zi), Z2 = 
diag(z2, Z2, Z2, Z2), pi(mi) = p{m)l4 and ^2(^2) = A(Z)l4. In particular, if Zi and 
Z2 are replaced respectively by e'^^UiZiUl and e^'^U2Z2U\ we have the results of 
Eq. 12.231 The special case p{m) = ml and A(/) = II may also be of interest. 

For an explicit demonstration, in the following example of the physical model 
we construct MVCS for a Hamiltonian, which is taken as the tensor product of 
two simple cases of the Jaynes-Cummings Hamiltonian (JC). As a second example, 
we construct VCS for the Hamiltonian (|1.8|) with the spectrum (|1.9|) shifted in a 
suitable manner. 

6.1. Tensored Jaynes-Cummings model. Let us construct MVCS for the Hamil- 
tonian H — ® He, where is the JC in the week coupling limit and He is the 
JC with the exact resonance and zero coupling (as described below). Information 
regarding the JC model is extracted from We outline it according to our needs. 
For a detailed analysis please see [T^ . 

The JC model describing a spin-i system in interaction with a one-mode magnetic 
field, in the rotating wave approximation, can be given by 

(6.5) Hjc = tj(a+a" + i)cro + ^0-3 + K(a+a-_ + a^<T+) 
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where uj is the field mode frequency, ujq the atomic frequency, k, a coupling constant, 
a+ and are the usual creation and annihilation operators for the radiation field, 
CT± and (73 are associated with the usual Pauli matrices and ctq is the identity matrix. 
The Hamiltonian Hyj: It is known that the Hamiltonian Hjc can be diagonalized 
as 

where O is the diagonalizing operator. From the diagonalized form the energy 
eigenvalues can be obtained as 

e~ = ujn + Kr{n) and = uj{n + 1) — Kr(n + 1) 

where r{n) = \/S + n, S = (^)^ and A = a; — wo, the detuning with A > 0. Since 
A > we have e^_^_i >e~. If < k/lu < 2^/6 + 1, the energies e+ are strictly 
increasing and non-degenerate. Let 

^±(k) = ^ and = e„ - €(, . 

In the weak coupling limit case we expand and by keeping at most terms of 
order 2 in k we get 

E^{k <<<) = LU±{K)n. 
Let be the corresponding normalized energy states. In this case we have 

p±{n)^E^Ef...Et = [L.±M]"r(n+l). 

The Hamiltonian He'- In the absence of the oscillating component of the magnetic 
field (k = 0) and for the exact resonance (A = — ojq = 0) we take Hjc '■= He- In 
this case it is straightforward to see that 

E^ = ijjn and p±{n) = w"r(n -f 1). 

Let (jy^ be the corresponding normalized energy states. 

For j = 1, 2, let z{,z{ e C and z{ r{e< , fc = 1, 2; j = 1, 2. Let 

Zi = diag(zi, Z2), Z2 = diag(zi, 

Pi{mi) = diag(p+(mi),p_(mi)), p2(?Ti2) = diag(p-^(m2), P-(to2)) 
and pi(0) = P2(0) = l2- Since, for k — 1,2, 



°° l'„l\2mi °° /^2\2m2 

y and y 

.1=0 \/P±("^l) „tto VP±("^2) 



converge everywhere we have — = cxd. From the inverse Mehin transform |18| 

poo 

Jo 



and Ht).2|l we observe that 



\\[r\) = 2rl exp ( A^Cri) = 2^^ exp f-^"^)^ 



A?(r^) = 2r^ exp ( _Ml!) and \l{rl) - 2r2 exp ^ ^''2)' 
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In this case a set of MVCS for H can be written as follows: 
(6.6) 



OO OO 



mi— m2—0 

k = 1,2, where {'&mi} and {$m2} are bases for two abstract separable Hilbert 
spaces Sji and Sj2, 

N{Zi,Z2) = cxp — — H H \ 

and a resolution of the identity is satisfied as 

2 f'oo POO POO POO p2tt p27v p2-n p2tv 



^ POO POO POO POO PZTT PZTT P Z7T P ZlT 

y, \Zi,Z2,k){Zi,Z2,k\dfi{ZuZ2) 
k=l "'0 "'0 "'0 "'0 "'0 Jo Jo 

= 12® Isji «) 

where 



2 2 



dM^i, ^2) = ——2-, n n ^'^'^^i- 



+ 3 = 1 k=l 

As a last illustrative example of the construction let us introduce VCS for the 
Hamiltonian H1.8|l with the spectrum (|1.9|l . 

6.2. The Hamiltonian (|1.8|) . Consider the Hamiltonian H of (|1.8I) with 5 = 
and denote it by Hq. In this case, the spectrum of p.9|) reads 

= + 1) + t^2(TO + 1) = e± + e,t 

Let ■i/'n ^ (/i^ be the corresponding normalized eigensolutions. Set 

R{n,m) = diag(£;!^^"!,i?!!^"M) = c^Ti^^Uh, 

where a = I + + 1) + LU2]/i-tJ2 and (a)„i = r(TO + a)/r(Q;) is the Pochhammer 
symbol. Let 

R{n) = diag(e+!,e-!) =r(n + 2)tjj'l2 
Z = diag(zi,Z2) and 
3 = diag(w,iJ) 

where zi = re'^i,Z2 = r2e^^'^,v = se"\v — se"''' and 61,92,1] € [0, 27r), where 
the matrix 3 is chosen in a convenient way (see below). As before, let {0n}^oj 
{'4'm}m=o orthonormal bases of two separable abstract Hilbert spaces Sji and 
f)2 respectively and {xi: X2} be the natural basis of C^. Identify the eigensolutions 
of -ffo to the basis of the abstract Hilbert spaces as follows: 

^^^^^ ^ f 



Q y := Xi «) 0« «) "0™ ^_ ^ ^_ j := X2 «> 0« ® -0™- 

With the above setup we write a set of VCS for the Hamiltonian Hq as follows: 

00 00 

I z, 3, k) = Miz, 3)-^ Ri^r^z^Mo, n)-^ J2 ^("' ^y^yxk ® 0™ ® 



n— m— 
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k — 1,2, where A/'(3, n) normalizes the second sum as introduced in (see Eq. 11.51 
or |K1 The normaUzation requirement 



k=l 



yields that 

oo 

where 



71=0 



w5T(n + 2) w5"r(n + 2) 



m— 



A/'(3,n)= ^ 



-^9 



iFi(l;a;— ) 
a;2 



converges everywhere. Since iFi(l; a; ^) > 1 for all s > 0, we have 



< 



E 



rt=0 



2n 



„2n 



Jin 



^T(n + 2) uj1T(n + 2) 

^2n 



^w^T(n + 2) w'j^T(n + 2) 



which converges everywhere. Let 

2r2 



o2q-1 



A(r) = — e- 



tU2"r(a) 



Then we have 



r^"A(r) / s^"A(s,n)dsdr = ^5^(71 + 2)^^(0)™. 
Jo 

Let D ~ [0,cx)) X [0, 27r) then by setting n — Q and m = in the above integrals 
one can see that 
1 



I DxDxDxD 

Thereby we have a resolution of the identity 
2 



A(ri)A(r2)A(s, n) drid0idr2d92dsdr/dsdri = 1. 



E 

k=l 



I Z, 3, k) (Z, 3, A: I d^liZ, 3) = I2 ® ® /i52 



DxDxDxD 



where 



1 

dAi(^,3) = j—^Y[jV{Z,5)X{rj)JV0,n)X{s,n)'^drjddjdsdsdT]d'n. 
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Here again Remark lfi.ll is applicable. For the Hamiltonian Hq we can have another 
set of CS similar to the CS (|6.6|l by appropriately shifting the spectrum. We 
describe it in the following remark. 



Remark 6.3. Let us shift the spectrum as follows: 

e^' = — = e J = -C/_ — hj = uj21ti 

Let 



n.m Tvn.m i-^O.n n.rn j^^n.ni T-iO.n 



R{n,m) = diag(e"'™!,e!!'"'!) =m!c<j™l2 
= diag(e!^!,e!l!) =n!<l2 
Z = diag (2:1,22) and 
3 = diag(wi,i;2) 

where zi = re*^i,Z2 = r2e**^wi = sie'''Sf2 = S2e*''^ and 0i,O2,r]i,r]2 € [0, 27r). 
With the above setup we write a set of VCS for the Hamiltonian Hq as follows: 

I Z,3,fc) =AA(Z,3)-^£^(")"^^" £ i?Km)-^3"Xfc0 0„<8)^™, = 1,2. 

n— m— 

By following the details of the Subsection IB . II we can obtain the normalization factor 
and the resolution of the identity. 



7. Conclusion 

The results of ^ have been generalized with multiple matrices. The results pre- 
sented in this article serve as a generalization of 'F since they can be employed for 
the purpose of obtaining CS for tensored solvable matrix Hamiltonians as well as to 
obtain CS for multi-level multi-mode solvable systems. Further the results of Sec- 
tion|Slmay also be considered as a generalization of the physical model presented in 
[7]. In order to demonstrate the physical model, in Example lG.ll we have considered 
the tensor product of two simple cases of the Jaynes-Cummings model. The model 
can also be applied to other tensored systems, where to solve the moment prob- 
lem one may adapt the results or the methods of Furthermore, the methods 
presented here may facilitate the derivation of CS for super Hamiltonians with mul- 
tiple degrees of freedom. The complex and real representation of quaternions and 
the real representation of octonions are used in the construction of multi-matrix 
VCS in several ways, namely the variables JjS of H1.5|l are replaced by matrices, the 
Pj{nj)8 are replaced by matrices, within the multiple sums one sum is allowed to 
depend on the other. In the same spirit, with the aid of 7 , it is possible to extend 
the construction to the matrix representation of any Clifford algebra and to any 
Cartesian product of classical matrix domains. 
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